Effects of magnetism and doping on the electron-phonon coupling in BaFe 2 As 2 
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We calculate the effect of local magnetic moments on the electron-phonon coupling in BaFe 2 As 2 +5 
using the density functional perturbation theory. We show that the magnetism enhances the total 
electron-phonon coupling by ~ 50%, up to A < 0.35, still not enough to explain the high critical 
temperature, but strong enough to have a non-negligible effect on superconductivity, for instance, 
by frustrating the coupling with spin fluctuations and inducing order parameter nodes. The en- 
hancement comes mostly from a renormalization of the electron-phonon matrix elements. We also 
investigate, in the rigid band approximation, the effect of doping, and find that A versus doping does 
not mirror the behavior of the density of states; while the latter decreases upon electron doping, the 
former does not, and even increases slightly. 

PACS numbers: 63-20.dk,63-20.Kd,74-20.Pq,74-70.Xa 



The simultaneous presence of high T c superconduc- 
tivity and magnetism in the phase diagram of the Fe- 
based superconductors (FBSC) suggests that magnetism 
plays an important role in determining the superconduct- 
ing properties. Phonons have been excluded early on as 
possible mediators, on the basis of first-principles non- 
magnetic calculations for the undoped compound, — ^ but 
the experimental situation is still far from settled in this 
regard. — ~— 

Density functional (DFT) calculations correctly repro- 
duce several properties of Fe pnictides, such as the mag- 
netic pattern in the parent compounds and the geometry 
of the Fermi surface. The interplay of magnetic and elas- 
tic properties is however puzzling— on one hand, exper- 
iments measure weak magnetic moments (m ~ 0.3 — 1.0 
Hb) in the SDW state and no long-range magnetic or- 
der in the superconducting samples, while spin-polarized 
LSDA calculations predict large magnetic moments at 
all dopings (m ~ 2.0 /is). On the other hand, the 
same spin-polarized LSDA calculations predict equilib- 
rium structures and phonon densities of states that are 
much closer to the experiment than those predicted by 
non-magnetic calculations. — ~— A possible way to recon- 
cile these two apparently conflicting results is an itinerant 
picture, in which the Fe atoms nevertheless have large lo- 
cal magnetic moments that order below the Neel temper- 
ature in the undoped compound, but survive locally at all 
dopings i 10 ' 11 . But this point of view immediately raises 
the question of whether the estimate of electron-phonon 
(e — ph) matrix elements given by early non-magnetic 
(rather than paramagnetic, i. e., with the local Fe mo- 
ment entirely suppressed) DFT calculations is represen- 
tative of the actual compounds. 12 

In this paper, we calculate from first principles the 
e—ph coupling constant in antiferromagnetic (AFM) and 
non- magnetic (NM) BaFe2As2, using the linear response 
method— £ We confirm that magnetism strongly affects 
the phonon frequencies, leading to a renormalization of 



the modes that involve Fe-As vibrations,— ~— and also find 
a strong effect on the e—ph matrix elements, leading to 
a ~ 50% increase with respect to the NM values. Using a 
rigid-band model, we show that the e—ph coupling con- 
stant A as a function of doping does not follow the density 
of states. Finally, we estimate the e—ph coupling in the 
paramagnetic (PM) state, by combining the nonmagnetic 
band structure (eigenenergies and eigenfunctions) with 
the magnetic phonon spectra and self-consistent poten- 
tials. For the relevant values of doping, we estimate an 
upper bound to the e—ph coupling constant A = 0.35, 
i.e. not high enough to explain superconductivity, but 
not sufficiently weak to be neglected. For instance, e—ph 
interaction may be one of the factors responsible for ex- 
perimentally observed gap anisotropy. 
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TABLE I: Calculated properties of BaFe2As2. m is the in- 
tegral over the cell of the absolute value of the magnetiza- 
tion, N(0) is the DOS in states/(spin eV Fe atom), \ aa 
is the electron-phonon coupling and uji„ is the logarithmic- 
averaged phonon frequency. In addition to the three fully 
self-consistent calculations we report here three model cal- 
culations (see text): PM1 utilizes the wave functions, one- 
electron energies and phonon frequencies from the NM calcu- 
lations, and the deformation potentials from the AFMc cal- 
culations. PM2 and PM3 use one-electron energies and wave- 
functions from the NM calculations, and phonon frequencies 
and deformation potentials from the AFMc and AFMs calcu- 
lations, respectively. 



2 



In order to disentangle the structural and mag- 
netic effects, we used the high-temperature tetragonal 
structur e 14 ! 15 in all our calculations. In order to es- 
timate the effect of the SDW ordering vector on the 
e — ph properties, along with the NM calculations, we 
have considered two different AFM patterns: the checker- 
board one (AFMc), in which the nearest neighbor Fe 
spins are antiparallel, and the experimentally observed 
one (AFM stripe, AFMs), in which spins are aligned 
(anti)ferromagnetically along the (y)x edge of the square 
Fe planes. The stabilization energies with respect to the 
NM solutions are 90 meV and 130 meV/Fe atom, respec- 
tively. 

All calculations were performed in the generalized gra- 
dient approximation 16 using plane-waves — and ultra- 
soft pseudopotentials 1 ^. We employed a cut-off of 40 
(480) Ryd for the wave- functions (charge densities). The 
electronic integration was performed over an 8 3 k-mesh 
with a 0.01 Ry Hermitian-Gaussian smearing in the NM 
and AFMc case, while for the AFMs case we used a 
8x4x8 k-mesh with a 0.01 Ryd Hermitian-Gaussian 
smearing. Finer grids (20 3 and 20 x 10 x 20 ) were 
used for evaluating the densities of states (DOS) and 
the phonon linewidths. Dynamical matrices and e — ph 
linewidths were calculated on 4 3 , NM and AFMc case, 
and 2 3 , AFMs case, uniform grids in q-space. Phonon 
frequencies throughout the Brillouin Zone were obtained 
by Fourier interpolation. The (perturbed) potentials and 
charge densities, as well as the phonon frequencies, were 
calculated self-consistently at zero doping (6 = 0); the ef- 
fect of doping on the e — ph coupling was then estimated 
using the rigid-band approximation. 

The band structures at <5 = (not shown) agree with 
previous calculations i 1 1 1 19 The corresponding magnetic 
moments and densities of states (DOS) at the Fermi level 
are reported in table HI 20 The phonon dispersion for NM, 
AFMc, and AFMs cases are shown in the left panels of 
Fig. [T] the remaining panels show the (partial) densi- 
ties of states and Eliashberg functions. Our results agree 
with previous calculations in the same crystallographic 
structures (adding orthorhombicity additionally changes 
phonon dispersions, see Refs. HH). Magnetism has the 
biggest effect on phonon modes that involve Fe-As vibra- 
tions; there is a pronounced softening of a branch origi- 
nally located at ~ 25 meV along the r — Z line in the NM 
calculation, down to ~ 20 meV in both AFMc and AFMs 
calculations. This branch corresponds to out-of-plane vi- 
brations of the As atoms. The in-plane Fe-As modes 
at low and high energy are sensitive not only to the size, 
but also to the pattern, of the magnetic moment. Indeed, 
when the order is AFMs, Fe vibrations along the AFM 
direction harden, whereas those along the FM direction 
soften, while for the As vibrations it is the opposite. 

The shift of phonon frequencies in the AFM calcula- 
tions has often been considered an indication of an en- 
hanced e — ph coupling in the magnetic phase ; 12 ' 24 but 
no explicit calculations of the e — ph coupling constant 



in the magnetic case have been reported so far. In this 
work, we have calculated from first principles the same- 




FIG. 1: (color online) Phonon properties of undoped 
BaFe2As2, for different magnetic patterns: non-magnetic 
(NM), AFM checkerboard (AFMc) and AFM stripe (AFMs). 
From left to right: the phonon dispersions; the total phonon 
densities of states, F(u}) (red, dashed), and Eliashberg func- 
tion, a 2 F(u>) (black, solid); partial contributions of Fe (red, 
dark) and As (green, light) ions to the total F(u)), projected 
onto Cartesian axes (the Ba contribution is limited to uj < 10 
meV, and not affected by magnetic order). The dispersions 
are shown in the same Brillouin zone for the three patterns. 
The k-points are selected so that they are physically the same 
in all three structures; the magnetic structure makes the two 
r — M directions inequivalent in the AFMs; the spins are 
aligned AFM(FM) in the x(y) direction. In this coordinate 
system, the points are: V = (0,0,0); Z = (0, 0, 7r/c) or 
(ir/a,ir/a,0); M = (tt/o, 0,0); M' = (0,7r/2a, 0), where a 
is the length of the shortest Fe-Fe bond and c is the Fe-Fe 
interlayer distance. 
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spin component of the Eliashberg spectral function: 
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Ck+qm = (k CT n|oV CT /<5e ql ,|k + q <7 m)/ v /2o; qiy . 



Here, iV^ is the number of k-points used in the summa- 
tion, N a (0) is the density of states per spin at the Fermi 
level, and Wq are the phonon frequencies. The e—ph ma- 
trix element g^nk+qm * s defined by the variation of the 
self-consistent crystal potential V a for the spin a with 
respect to a frozen phonon displacement according to 
the phonon eigenvector e q!/ = YLao. m a V^V 6 ^^ • 
Here UqAa is the Fourier transform of the a component of 
the phonon displacement of the atom A in the unit cell, 
Ma is the mass of atom A and e\ a are Aa components 
of qis phonon eigenvector normalized in the unit cell. 
The first inverse moment of a 2 F(u>) gives the frequency- 
dependent e—ph coupling constant: 



\ aa (u) = 2 / dna 2 F aa {n)/n. 



(2) 



The total e — ph coupling constant \ aa — A CTCT (w = oo) 
is 0.18,0.33 and 0.18 for NM, AFMc and AFMs, respec- 
tively. 

The three calculations in Fig. [T]have different phonon 
spectra w^, different self-consistent crystal potentials 
V a , and different one-electron wave- functions |kn) and 
eigenenergies e£, which all determine the value of A in 
Eqs. [TJH In order to understand which of these factors 
dominates the e—ph coupling, we start for simplicity 
from the so-called Hopfield formula for the e—ph cou- 
pling constant: 



Xn 



N a (Q)D 2 
Moj 2 



(3) 



where D is the deformation potential, which is a measure 
of the average e—ph interaction, and Mur is a character- 
istic force constant, which we assume depends weakly on 
the magnetic order. We can then use the ratio X aa /N a (0) 
to obtain an estimate for the average change in the e—ph 
interaction due to magnetism; going from NM to AFMc 
(AFMs), this ratio, and hence D 2 , increases by ~ 50%, 
from 0.15 to 0.24 (0.26) - see table Q] A comparison of 
Eq. [T] and Eq. [3] shows that the increase in D can be 
caused either by an increase in the derivative of the DFT 
potential V a with respect to displacement, or by a differ- 
ence in the one-electron wavefunctions, used to evaluate 
its average over the Fermi surface. In order to disentan- 
gle these two effects, we performed a mixed calculation 
(PM1), in which we combined the NM eigenvalues and 
wavefunctions, as well as the phonon frequencies, with 
the AFMc potential variations in the Eq. [TJ In the top 
panel of Fig.[2J we compare the resulting Eliashberg func- 
tion a 2 F(td) with the NM one; the total phonon DOS 
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FIG. 2: (color online :) The Eliashberg functions a 2 F(uj) 
for the non- magnetic (NM, top, black solid line), anti- 
ferromagnetic (AFM, middle and bottom, black solid lines) 
and model paramagnetic calculations, as described in the text, 
(PM1-3, colored dashed lines). The corresponding phonon 
densities of states, in arbitrary units, are shown in each panel 
by the dotted lines. 



F(oj) is also shown as dotted lines. The effect of us- 
ing the AFM potential is an increased coupling of elec- 
trons to phonons with frequencies u ~ 20 — 25 meV, 
while other modes are largely unaffected. This increases 
A from 0.18 (NM) to 0.27 (PM1). The middle and lower 
panel of Fig. [5] show two calculations, PM2 and PM3, 
which use NM one-electron wave functions, and AFMc 
and AFMs phonon frequencies and crystal potentials, re- 
spectively. A comparison of the a 2 F(ui) with the phonon 
densities of states shows that, also in these cases, there is 
an increased coupling to phonons around 20 meV, both in 
AFM and PM. The A's, reported in table H are 0.27 and 
0.31, for PM2 and PM3 respectively. From the compari- 
son of PM1 and PM2, we conclude that the effect of the 
phonon frequency is negligible, while the ~ 10% spread 
of values between PM2 and PM3 gives an indication of 
the effect of the short range AFM correlations on A. Fi- 
nally, we find that the values of X aa /N a (Q) in the PM 
calculations are in line with the AFM ones, indicating 
that, at S = Q, the states at Ep have a comparable, weak 
coupling to phonons, for both NM and AFM orders. 

These model calculations, which combine the AFM 
potentials and phonons and the NM wavefunctions and 
Fermi surfaces, represent the best approximation, at 
LDA level, of the real paramagnetic (PM) state of su- 
perconducting, doped, BaFe2As2 which is characterized 
by local, disordered magnetic moments^ In the follow- 
ing, we will use them to estimate an upper bound for the 
e — ph coupling in the doped BaFe2As2- 
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FIG. 3: (color online) Effect of doping on the e—ph coupling 
of BaFe2As2, in the rigid-band model. ±<5 is the number of 
excess electrons/holes per Fe atom. The left panel shows the 
density of states at Ef, in states/spin eV Fe atom. The right 
panel shows the same-spin coupling constant A = A CTCT . N a (0) 
is the same for NM and all three (see text) PM calculations, 
as described in text. 

In Fig. wc show the density of states af the Fermi 
level N a (0) (left) and \ aa (right) as a function of doping 5 
for NM, AFMc, AFMs, and PM orders, from a rigid-band 
calculation; i5 is defined as the number of excess electrons 
(holes) per Fe atom. NM and PM calculations have the 
same -/V^O), which decreases monotonically as a function 
of doping, while \ aa is roughly symmetric around 5 = 0. 



This indicates that, in Eqs. [T][3l the effect of the matrix 
element dominates over that of the DOS. In the AFM cal- 
culations, the shape of the Ao-o- curve follows more closely 
that of N a (0), with a minimum at 6 ~ 0.1 — 0.2, and a 
maximum \ aa = 0.8 for AFMs at 8 = —0.5, correspond- 
ing to KFe2As2- However, the use of the rigid-band ap- 
proximation is questionable at these high dopings, and in 
the following we limit our analysis to a smaller range of 
5 = ±0.25, where the error in A connected to the rigid- 
band approximation is < 20%^ 

We can summarize the results of Fig. [3J by saying 
that, for 1 5 < 0.25, the upper bound for the e — ph 
coupling in paramagnetic, superconducting BaFe2As2 + 
2<5 is A (j (j < 0.35. The enhancement in A results from a 
~ 50 % increase in the e — ph matrix elements, due to 
static magnetism, which is independent on doping, and 
a symmetric increase of the matrix element for 6 =/= 0, 
which does not follow the shape of the electronic DOS. 

The value X aa = 0.35, which includes magnetism and 
doping, is almost a factor of two larger than that esti- 
mated in early non-magnetic calculations for the undoped 
compounds.ii 2 >, and close to to a recent experimental es- 
timate of A from a kink in photoemission spectral 

What are the consequences of this result? Of course, 
it cannot explain a T c of 38 K, which confirms that su- 
perconductivity in this and other Fe-based superconduc- 
tors is most likely due to electronic (magnetic) degrees 
of freedomi2i2£ However, in Fe pnictides, studies of the 
superconducting gap in realistic models with AFM spin 
fluctuations show that solutions with and without gap 
nodes are almost degenerate, so that even a relatively 
low e — ph coupling constant can help select either of 
them, by enhancing/suppressing the pairing in the rele- 
vant channeL 24 ' 25 
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